Loss separation for dynamic hysteresis in magnetic thin films 
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We develop a theory for dynamic hysteresis in ferromagnetic thin films, on the basis of the 
phenomenological principle of loss separation. We observe that, remarkably, the theory of loss sep- 
aration, originally derived for bulk metallic materials, is applicable to disordered magnetic systems 
under fairly general conditions regardless of the particular damping mechanism. We confirm our 
theory both by numerical simulations of a driven random-field Ising model, and by re-examining 
several experimental data reported in the literature on dynamic hysteresis in thin films. All the 
experiments examined and the simulations find a natural interpretation in terms of loss separation. 
The power losses dependence on the driving field rate predicted by our theory fits satisfactorily all 
the data in the entire frequency range, thus reconciling the apparent lack of universality observed 
in different materials. 
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PACS numbers: 76.60.-d,75.60.Ej,75.70.-i,75.70.Kw 

Power losses in ferromagnetic materials generally de- 
pend on the frequency of the applied field, a phenomenon 
referred to as dynamic hysteresis 0, 0]. The problem 
has important implications for applications to high fre- 
quency devices, and, from a purely theoretical point of 
view, for the understanding of the dynamics of disor- 
dered magnetic systems, which represents a central issue 
in non-equilibrium statistical mechanics. While dynamic 
hysteresis in metallic bulk three dimensional systems is 
well understood in terms of eddy current dissipation, a 
satisfying theory for thin films, where the effect of eddy 
current is expected to become negligible 0], is still lack- 
ing. Hence^in recent years a g reat attention, both exper- 
imental [lIlllEMMEflElElElElElE^ and 

theoretical [l|, [lT], Il8l [l9j , has been devoted to magnetic 
reversal in thin and ultrathin ferromagnetic films. 

An accurate interpretation of the experimental data 
requires a detailed understanding of the magnetization 
reversal properties on a microscopic scale. Based on the 
analogy with Ising type models, experimental data are 
often analyzed in terms of universal scaling laws, such as 
the one relating the dynamical hysteresis loop area (A) 
to external parameters such temperature (T) , amplitude 
(Hq), and frequency (to) of the applied magnetic field: 
A oc HqUj^T^ 1 , where a, 0, and 7 are scaling exponents 
0. The experimental estimates of these exponents, of- 
ten based on a very limited sca ling reg ime, display a huge 
variability 0,0,0,0,0,0,113,11111131, so that the valid- 
ity of a simple universal scaling law is still under ques- 
tion. Some authors also interpret the lack of good scaling 
of A{io) as a cross-over between two distinct dynamical 
regimes, at low frequencies dominated by domain wall 
propagation, at high frequencies by the nucleation of new 
domains @,E3 On the other hands, several authors rec- 
ognized the present of a static hysteretic term, yielding a 
non-vanishing loop area in the limit w — ► [la, llfil |2fJ • 
While this surely represents an important step forward 
for a better description of the data, we still need a more 



comprehensive theory to understand the magnetization 
reversal in thin films. 

For bulk ferromagnetic materials, the conventional the- 
ory of power losses takes as a starting point the phe- 
nomenological principle of loss separation according to 
which the average power loss can be decomposed into 
the sum of static (hysteretic), classical and excess contri- 
butions: P = Physt +Pci + Pex 01 The hysteretic contri- 
bution Physt is due to the presence of quenched disorder, 
and gives rise to a non vanishing area of the quasi-static 
hysteresis loops. The classical term P c i is obtained com- 
puting the eddy currents present in an equivalent homo- 
geneously magnetized sample j^. Finally, excess losses 
P ex take into account other contributions due to the dis- 
sipation associated to the presence of multiple moving 
domain walls • The theory of loss separation is usu- 
ally considered inadequate for thin films, since eddy cur- 
rents dissipation disappears as the sample thickness goes 
to zero. While this is certainly correct for what concerns 
the contribution from classical losses, the hysteretic and 
excess loss contributions have been derived under very 
general conditions which are in principle valid also for 
thin films Q. 

In this letter, we show that once the excess loss contri- 
bution is taken into account, the theory of loss separation 
describes very accurately the properties of dynamic hys- 
teresis in thin films. Instead of directly considering power 
losses, we focus on the behavior of the coercive field H c 
as a function of the driving field rate H = dH/dt, since 
these are the quantities used in published experimental 
data. Analyzing the behavior of coercive field is particu- 
larly convenient also because it is less sensitive than the 
loop area to the particular driving condition employed 
in the experiment (i.e. sinusoidal or triangular field, ap- 
plied field or flux rate, etc.). Moreover for fully saturated 
loops the two quantities are proportional. 

Analogously to the total losses, the coercive field can 



be separated into three contributions 
H c = H„ + H c j + H ea 



(1) 



where H p is the static (hysteretic) contribution of the 
pinning field, H c i is the classical eddy current contribu- 
tion, and H ex is the excess field. 

The classical term H c \ — C C \H is linear in the driving 
field rate H through a coefficient C c / oc cPafj, 3J, where 
a is the conductivity, fi is the permeability, and d is the 
sample thickness. Clearly, this term goes rapidly to zero 
with d, and thus it can be safely disregarded in thin films 
and in the following we set C c i = 0. 

The excess field has no explicit dependence on the 
thickness, thus in principle it should be taken into ac- 
count also in thin films. Originally, it has been derived 
under two basic hypothesis: i) eddy currents provide the 
dominant dissipation mechanism of magnetization rever- 
sal, and ii) a number of correlated regions between mag- 
netic objects (domain walls, for instance) are progres- 
sively activated by the increase of the the magnetization 
rate /. As a matter of fact, these hypothesis can be 
made more general. If we consider a single magnetic ob- 
ject which could reverse all the sample magnetization, the 
total excess field H w is taken proportional to the magne- 
tization rate /, so that H w = VI, where T is an effective 
damping coefficient, This coefficient was originally calcu- 
lated for the eddy current dissipation, but that can be 
generalized to any general dissipation mechanism (spin 
relaxation, etc.) which produces an excess field H w . Be- 
cause in thin films the control parameter is H, we also 
assume / = /j,H, where fj, is the permeability. Given that, 
the factor T/i turns out to be an effective relaxation time 
relating the applied field rate to the excess field for a 
single magnetic object, so that H w = YfiH. 

In the case of n magnetic objects, the excess field H ex 
is proportionally reduced giving H ex — H w /n. At the 
same time, the number of active magnetic objects itself is 
a function of the excess field, and, as in the original model 
|22j . we assume a simple linear relation n ~ no+H exc /Vo, 
where no is the number of active magnetic objects in the 
quasi-static limit, and Vb is a characteristic field which 
controls the increase of the number of active magnetic 
objects due to the excess field. 

With these hypothesis, the excess field is calculated as 



C ex ((l + rH)^ 2 -1) 



(2) 



where C ex = n Vo/2, and r = 4T(j,/(iIqVq). The excess 
field thus only depends on the typical relaxation time 
Tfi, the number of magnetic objects in the quasi-static 
limit no, and the characteristic field Vq. Remarkably, as 
reported in when rH <C 1, we get H ex = T^,H/n , 
that is, for a large number of magnetic objects (n ~ no 3> 
1) the excess field is a linear function of the applied field 
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FIG. 1: Variation of dynamic coercive field H c with the ap- 
plied field rate dH/dt, for a Fe-based amorphous ribbon with 
thickness 20 fim. The coercive field measured under constant 
field rate H (squares) and constant magnetization rate / (cir- 
cles) collapses using H = I / (J,, where u = 0.08 is the linear 
permeability around the coercive field. 



i.e. the excess field has a square-root dependence on the 
rate when the number of magnetic objects is pretty small 
(no~l). 

At this point, one can argue about the equivalence of 
performing experiments by controlling the driving field 
H instead of the magnetization rate / in dynamic exper- 
iments. To verify this assumption, we have performed dy- 
namic hysteresis experiments in an Fe^Co^H^ amor- 
phous ribbon having a thickness of 20 \xm, a typical sam- 
ple where the loss separation is known to apply. We have 
measured the coercive field both as a function of H and 
of /. As shown in Fig. ^ the data from the two exper- 
iments collapse reasonably well using the linear perme- 
ability /i ~ 0.08 measured around the coercive field. 

In order to test the generality of the loss separa- 
tion mechanisms for hysteresis we consider the driven 
random-field Ising model (RFIM), which has been pro- 
posed in the past as a paradigmatic model to under- 
stand the effect of disorder in ferromagnetic hysteresis 
El EH In the RFIM, a spin s l = ±1 is assigned 
to each site of two dimensional square lattice. The spins 
are coupled to their nearest-neighbors spins by a ferro- 
magnetic interaction of strength J and to the external 
field H . In addition, to each site is associated a ran- 
dom field hi taken from a Gaussian probability density 
p(h) — exp(— h 2 /2R 2 )/y/2irR, with width R, which mea- 
sures the strength of the disorder. The Hamiltonian thus 
reads 



n 



E 



(H + hi)si 



(3) 



rate. On the other hand, for rH 3> 1, H e 



TfxV H, 



where Yliij) ^ s restricted to nearest-neighbors pairs. We 
consider a simple relaxation dynamics obtained in the 
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FIG. 2: The dynamic coercive field as a function of the applied 
field rate from simulations of the driven RFIM with different 
values of the variance R of the disorder distribution. The 
curves are fitted according to Eq. |2] with C c ; = 0. In the 
inset: data in lin-lin scale. Note the linear behavior at large 
R, and the square-root dependence at low R. 



limit T -> of the Glauber dynamics [23] |2J, |2g : at 

each time step the spins align with the local effective 
field Si = sign ^ J £\ Sj +h t + Hj. 

This model was originally studied in the quasistatic 
limit, however, dynamic effects can easily be incorpo- 
rated assigning a timescale r to spin relaxation and in- 
creasing the magnetic field in steps A. This is equivalent 
to an effective field rate H = A/r. We numerically com- 
pute hysteresis loops for different values of the disorder 
strength R and obtain the coercive field averaging over 
different realizations of the disorder. In Fig. |3 we show 
that the data are described well by EqsnjJ2l The result- 
ing fit indicates that C ex grows with the disorder, while 
r decreases. This fact can be understood by noting that 
the growth of disorder naturally leads to an increase in 
the number of the domains. In fact, when R ~ 1, no is 
pretty small so that rH ^> 1 and the excess field has a 
square-root dependence on the rate H. At R ~ 3, no 
is large (rH < 1) and the dependence is linear (see the 
inset of Fig.EJ. 

It is interesting to remark that in Ref. |2fJ similar re- 
sults for the loop area in the RFIM have been fitted by 
a scaling law A = Aq + A\lo@ ', obtaining a disorder de- 
pendent (3. While it is difficult to discriminate between 
different fitting functions, we notice that an exponent de- 
pending on microscopic parameters is difficult to justify 
theoretically. In this respect, Eq. ^ interpolates between 
(3 = 1 for rH < 1 and f3 — 1/2 for rH > 1, as said. We 
also note that (3=1/2 was obtained from the solution of 
a model for the dynamics of a single domain wall |l6ll27| . 
fully compatible with the present results. 

As Eqs. provide a good description of dynamic 
hysteresis under fairly general conditions, we can try to 
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FIG. 3: Analysis of experimental data using Eqs. [TH^lon epi- 
taxial Fe/GaAs(001) thin films, as reported in Ref. Data 
for the 25 nm film (■) are approximately described by a linear 
function 
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FIG. 4: Analysis of experimental data using Eqs. EE1 on a 
MoS2/Au/Co/Au sandwich taken at 300 K, as reported in 
Ref. [3. The excess field has approximately a square- root 
dependence on the applied field rate. 



compare its predictions with some of the experimental 
measurements reported in the literature. In particular, 
we re-examine the measurements of dynamic coercivity 
in magnetic thin films and multilayers of various thick- 
nesses, as reported in Figs. (from Refs. 0;E0,E|], 
respectively). Significantly, all the data are well fitted 
by Eqs. over the entire range of the applied field 
rate. Both limiting cases are present (rH <C 1, and 
rH 3> 1), with an approximately linear (black squares) 
and a square-root (black triangles) dependence on H. 
Despite the large variability of compositions, structures, 
and static coercive fields (which span two orders of mag- 
nitude), the fitting curves are compatible with a lim- 
ited range of the model parameters: no = 1 — 100, 
Vb = 0.1 - 10 Oe, and Tfi = 10~ 2 - 10~ 4 s. In particular, 
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the time costant lays within the expected range for the 
dynamics of domain walls in thin films 0] . This observa- 
tion strongly supports the validity and the large applica- 
bility of this model. This also seriously puts in doubt the 
existence of a dynamic transition between a regime dom- 
inated by propagation and one dominated by nucleation 
of new domains. As a matter of fact, the model does not 
exclude the nucleation of domains, as an increasing value 
of n can imply both the activation of existing magnetic 
objects and the creations of new ones. But we also must 
note that in the case of a square-root dependence on the 
rate, the number of magnetic objects cannot change sig- 
nificantly, and that also in the opposite case of a linear 
dependence, when the number of static magnetic objects 
is large the effect of nucleation of new active domains is 
practically irrelevant. 
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FIG. 5: Analysis of experimental data using Eqs. on 
a NiFe and Co single magnetic layer films, as reported in 
Ref. fl^l . Data for the NiFe (■) are approximately described 
by a linear function. 

In summary we show that the properties of dynamic 
hysteresis in ferromagnetic thin films can be explained in 
terms of the theory of loss separation, originally derived 
for bulk metallic materials. This conclusion is supported 
by numerical simulations and by the analysis of experi- 
mental data reported in the literature. 
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